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ABSTRACT 


The  use  of  boiling  heat  transfer  in  the  liquid  immersion 
cooling  of  electronic  components  has  always  been  hampered  by 
the  excessive  superheat  necessary  to  initiate  nucleation  in 
the  fluorinated  hydrocarbons  used  as  dielectric  cooling 
fluids.  In  an  attempt  to  overcome  some  of  these  difficulties, 
an  experimental  study  of  the  effects  of  nucleate  pooling 
boiling,  on  the  heat  transfer  from  surface  near  the  boiling 
surface  was  conducted.  An  experimental  chamber  was  constructed 
which  had  a  column  of  four  horizontal  wires  spaced  2.5  cm 
vertically  from  each  other.  The  lowest  wire  was  progressively 
heated  from  the  natural  convection  region  through  the  nucleate 
boiling  region.  A  study  was  made  of  the  effects  of  the  boiling 
wake  from  the  lowest  wire  on  heat  transfer  from  the  upper 
wires.  Under  certain  conditions  heat  transfer  enhancements  of 
up  to  30%  were  obtained. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Research  on  pool  boiling  in  electronic  cooling  systems  has 
focused  on  three  primary  areas:  (1)  reducing  the  temperature 
excursion  at  incipient  boiling,  (2)  reducing  wall  superheat 
during  nucleate  boiling  and  (3)  enhancing  critical  heat  flux 
(Mudawar,  1990) .  Since  the  heat  transfer  potential  of  nucleate 
boiling  is  well  known,  this  regime  was  one  natural  selection 
for  examination  as  a  possible  solution  to  this  cooling 
problem.  Nishikawa  (1980)  notes  that  there  are  two  primary 
methods  used  to  promote  nucleate  boiling.  The  first  is  to 
treat  the  surface  in  a  manner  that  reduces  its  wettability. 
The  second,  more  promising  one,  is  to  manufacture  a  surface 
with  cavities  which  trap  and  hold  vapor  keeping  the  nucleation 
site  active. 

In  recent  years  the  drive  toward  higher  heat  flux  and 
power  density  in  electronic  equipment  has  turned  attention  to 
the  use  of  dielectric  liquias.  The  3M  Corporation  has 
commercially  produced  its  "Flourinert"  series  of  high 
dielectric,  inert  electronic  cooling  liquids.  In  contrast  to 
water  and  other  conventional  coolants  these  liquids  possess 
relatively  low  thermal  conductivities  and  boiling  points  and 
small  heats  of  vaporization.  They  demonstrate  high  wettability 
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on  most  surfaces.  Their  surface  tensions  are  the  lowest  of  any 
known  liquids,  resulting  in  near-zero  wetting  angles.  This  is 
significant  because  it  eliminates  many  potential  nucleation 
sites.  High  dielectric  strength  and  high  resistivity  are  some 
of  the  characteristics  that  make  these  liquids  well  suited  for 
electronic  applications.  Especially,  FC-72  is  known  as 
chemically  stable  and  leaves  no  residue  during  boiling. 

However,  when  using  dielectric  fluids,  incipience 
superheat  excursion  is  pronounced  which  leads  to  very  high 
wall  temperatures.  A  firm  understanding  of  boiling  incipience 
is  required  in  order  to  implement  effective  cooling.  Boiling 
can  produce  very  large-scale  heat  transfer  coefficients  but  it 
also  introduces  some  problems:  (i)  Boiling  restricts  the" 
physical  design  of  the  system,  (ii)  A  high  degree  of  superheat 
may  be  required  if  the  surface  is  very  smooth  in  order  for 
nucleation  to  occur.  Due  to  the  high  local  heat  flux  the 
occurrence  of  periodic  boiling  may  result  in  local  turbulence, 
even  damage  of  electronic  components.  (iii)  Impurities 
concentration  may  result  at  the  nucleation  site. 

Subcooling  offers  many  advantages  over  saturated  boiling 
for  cooling  large  arrays  of  computer  chips  due  to  the 
decreased  diameter  and  rapid  collapse  of  departing  bubbles. 

In  this  study  heat  transfer  enhancement  due  to  bubble 
pumping  in  FC-72  was  investigated.  It  is  hoped  that  they  will 
serve  to  identify  some  trends  and  provide  a  basis  for  further 
research  on  actual  components  and  integrated  circuits. 
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B.  PREVIOUS  WORK 


Hsu  (1962)  proposed  a  method  for  determining  the  range  of 
active  nucleation  sites  as  a  function  of  temperature  or  heat 
flux.  This  model  revealed  that  the  maximum  and  minimum  sizes 
of  active  cavities  are  functions  of  subcooling,  system 
pressure,  physical  properties  and  the  thickness  of  the 
superheated  liquid  layer. 

Bergles  and  Chu  (1979)  compared  the  nucleate  pool  boiling 
heat  transfer  characteristics  of  three  copper  Union  Carbide 
High  Flux  test  sections  to  a  plain  copper  tube  in  distilled 
water  and  R-113.  The  experiments  were  conducted  with  three 
different  treatments  to  the  test  surface  prior  to  collecting 
data:  subcooling;  aging  by  preboiling  in  the  pool;  and  heating 
the  surface  in  air  to  remove  all  liquid  in  the  pores.  This 
research  indicated: 

1.  Heat  transfer  coefficients  for  the  porous  surfaces 
were  four  to  ten  times  higher  than  for  the  plain 
surfaces . 

2.  There  was  a  significant  temperature  overshoot 
prior  to  the  initiation  of  boiling  and  it  was  not 
very  sensitive  to  aging  or  the  power  increment 
changes  in  water.  These  overshoots  were  sensitive 
to  aging,  initial  subcooling  and  power  increment 
changes  in  R-113. 

Bergles  et  al.  (1968),  conducted  nucleate  pool  boiling 
studies  in  water,  R-113  and  FC-78,  using  a  304  stainless  steel 
tube.  As  heat  fluxes  were  increased,  the  heat  transfer 
coefficient  followed  the  convective  curve  until  a  high 
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superheat  was  attained,  and  then  a  distinct  increase  in  the 
heat  transfer  coefficient  was  noted  as  nucleation  occurred. 
This  was  followed  by  an  abrupt  drop  in  all  temperatures  by  as 
much  as  27°F.  It  was  also  noted  that  the  inception  of 
nucleation  could  be  triggered  by  vibrating  the  test  surface. 
The  conclusions  drawn  were: 

The  phenomenon  of  temperature  overshoot  hysteresis  in 
ordinary  liquids  is  due  to  two  causes: 

1.  The  existence  of  metastable  bubbles  which  are 
triggered  only  at  sufficiently  high  disturbance 
levels, 

2.  The  deactivation  of  larger  cavities  by  displacing 
the  vapor  by  liquid  during  subcooling. 

Many  of  the  experimental  reports  published  prior  to  1986 
are  reviewed  by  Bar-Cohen  and  Simon  (1988)  .  They  noted 
significant  variations  in  the  incipience  superheat  from  one 
report  to  the  next.  The  authors  discussed  possible  mechanisms 
for  delayed  nucleation.  Recently  non-repeatable  and  unsteady 
behavior  of  boiling  incipience  for  a  highly-wetting  dielectric 
liquid  (R-113)  was  documented  by  You  et  al.  (1989). 

The  bubble  behavior  from  the  nucleate  boiling  regime  to 
the  film  boiling  regime  has  not  been  clarified  yet,  except  for 
liquid  helium.  The  transient  boiling  experiment  of  liquid 
nitrogen  at  stepwise  heat  generation  in  a  platinum  wire  was 
carried  out  by  Kunito  Okuyama  (1989) . 
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McAdams  et  al.  (1949),  reported  the  effect  of  dissolved 
air  on  forced-convection  subcooled  nucleate  boiling  of  water 
where  the  dissolved  air  reduced  the  heat  flux  at  boiling 
initiation . 

Torikai  et  al.  (1970),  used  platinum  wires  immersed  in 
water  to  measure  dissolved  gas  effect  on  the  incipience  wall 
superheat  under  reduced  pressure  conditions  (0.05  -  1.0  Bar). 

Surface  enhancement  is  a  very  popular  tool  for  improving 
pool  boiling  performance.  Marto  and  Lepere  (1982)  found  three 
commercially  available  enhanced  boiling  surfaces  which  shifted 
the  boiling  curve  of  FC-72  toward  lower  superheats  compared  to 
a  plain  copper  surface.  In  a  study  by  Berenson  (1962), 
nucleate  boiling  heat  transfer  coefficients  varied  by  600%  due 
to  variation  in  surface  finish.  He  emphasized  the  importance 
of  the  surface  roughness.  Recent  investigations  on  surface 
effects  in  pool  boiling  by  Chowdhury  and  Winterton  (1985)  show 
that  heat  transfer  appears  to  be  unaffected  by  contact  angle. 

The  effects  of  flow  velocity  and  subcooling  on  critical 
heat  flux  from  a  simulated  electronic  chip  attached  to  the 
wall  of  a  vertical  rectangular  channel  was  determined  by 
Mudawar  and  Madox  (1988) .  They  concluded  that  heat  flux  was 
significantly  enhanced  by  increasing  the  degree  of  subcooiing. 

The  solubility  of  gas  in  highly-wetting  dielectric  fluids 
is  of  special  concern.  These  fluids  can  contain  up  to  25  times 
more  dissolved  air  by  volume  than  can  water:  48%  air  by  volume 
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can  be  dissolved  in  FC-72  compared  to  1.9%  air  by  volume  in 
water  (3M  Product  Manual,  1987) . 

Experimental  results  presented  by  You,  Simon  and  Cohen 
(1990)  indicate  that  dissolved  air  plays  an  important  role  at 
incipience  only  when  a  very  high  dissolved  air  content  is 
maintained.  The  surface  tension  value  for  FC-72  was  estimated 
to  have  a  2%  reduction  due  to  high  dissolved  air  content. 

S.M.  You  et  al.  (1989)  conducted  an  experiment  in 
saturated  R-113  fluid  under  atmospheric  pressure  with  two 
types  of  wires.  One  is  0.13  mm  diameter  chromel  wire  (90% 
nickel)  and  the  second  is  a  thin  platinum  film  on  a  quartz 
rod.  They  examined  many  pits  of  various  sizes  using  a  scanning 
electron  microscope  (SEM) . 

Bar-Cohen  (1990)  conducted  experiments  providing  pool 
boiling  characteristics  for  a  number  of  commercially  available 
inert  liquids  (FC-43,  FC-72,  FC-75,  FC-77,  FC-84,  FC-87)  .  Data 
were  obtained  with  a  classical  0.25  mm  platinum  wire  heating 
element.  He  concluded  that: 

1.  Nucleate  boiling  characteristics  of  these  fluids 
depend  on  the  heating  surface, 

2.  Saturated  pool  boiling  curves  were  almost 
identical  with  one  another  on  a  particular  heater, 

3.  Superheat  excursions  at  boiling  incipience  appear 
to  always  occur,  but  their  extent  can  vary  widely. 
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C .  OBJECTIVES 


The  objective  of  this  work  was  to  investigate  the  effects 
of  the  boiling  wake  from  a  heated  wire  on  the  heat  transfer 
from  wires  placed  within  the  wake.  This  work  was  carried  out 
using  the  dielectric  liquid  FC-72. 
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II.  EXPERIMENTAL  DESIGN 


A.  DESCRIPTION  OF  COMPONENTS 

Figure  (1)  shows  the  system  configuration  and  identifies 
all  major  components.  Figure  (2)  is  a  photograph  of  the 
apparatus . 

The  chamber  consisted  of  a  plexiglass  box  with  an  aluminum 
coverplate  fitted  with  a  rubber  0-ring.  This  vessel  has  side 
wall  heaters  to  allow  the  liquid  to  be  preboiled  for  degassing 
before  the  run,  and  to  maintain  the  liquid  at  saturation 
temperature  throughout  the  run.  The  vapor  was  condensed  and 
returned  to  the  chamber  by  gravity  from  the  aluminum  cover 
plate  condenser  equipped  with  thermoelectric  coolers.  The 
following  is  a  description  of  the  test  apparatus. 

1 .  Tost  Chamber 

The  test  chamber  was  constructed  of  plexi-glass  of  1/2 
inch  thickness  and  had  inside  dimensions  of  6  inches  height, 
6  inches  length  and  2  inches  width.  An  isometric  view  of  the 
chamber  is  presented  in  Figure  (3)  .  Side  walls  are  glued  to 
each  other. 

It  has  a  1  inch  diameter  hole  and  1  inch  protruding 
pipe  and  a  copper  elbow  tube  in  the  back  face  which  is  the 
pathway  for  the  power  and  thermocouple  wires.  The  chamber  also 
allows  the  test  section  to  slide  in  and  out  for  the 
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experiment.  Cross-sectional  and  schematic  drawings  of  the  test 
section  in  the  chamber  are  presented  in  Figures  (4)  and  (5) 
respectively.  The  chamber  also  contained  two  small  (1/4  in) 
pressure  equalization  outlets,  located  at  the  top-right  and 
the  bottom-left  of  the  front  vertical  surface.  It  is  a  simply 
a  closed  tank  with  an  aluminum  cover  plate. 

The  test  chamber  is  attached  to  a  wood  base  in  order 
to  give  better  arrangement  for  cables  and  connections. 

2 .  Thermocouples 

Thermocouples,  0.01  in  (0.254  mm)  in  diameter,  were 
used  to  take  measurements.  The  thermocouples  for  the  plate 
condensers  (thermoelectric  coolers  which  are  located  on  top  of 
the  aluminum  plate)  were  bonded  in  grooves  by  placing  a  small 
drop  of  Eastman  910  adhesive  on  the  thermocouple  bead.  The  two 
thermocouples  then  were  placed  in  the  designated  line  groove. 
Using  a  weight,  pressure  was  applied  until  the  adhesive  set. 
After  allowing  for  three  hours  cure  time,  the  remainder  of  the 
groove  was  filled  with  Omega  Bond  101  Epoxy  and  smoothed  to 
the  level  of  the  aluminum  plate  inside  face.  The  Omega  Bond 
101  Epoxy  experienced  little  shrink  after  curing. 

3 .  Heaters 

The  side  wall  heaters  were  glued  on  either  side  walls 
of  the  chamber.  The  actual  size  of  the  heaters  were  125  x  11 
mm  and  they  were  subjected  to  5  V,  2.5  amperes.  The  power  load 
depended  on  the  combination  of  the  wires  being  heated. 
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4 .  Boat  Exchangers 

The  three  thermoelectric  heat  exchangers  measure  1.56 
X  1.56  inches  and  were  placed  on  top  of  the  aluminum  cover 
plate.  The  power  supply  for  these  were  determined  by  the 
related  curves.  The  property  curves  and  technical  data  are 
presented  in  Figures  (6)  and  (7)  respectively. 

5.  Aluminum  Cover  Plata 

An  aluminum  plate  of  1/2  inch  thickness  was  used  as  a 
cover  plate  for  the  test  chamber.  The  inside  face  of  the  plate 
had  two  thermocouple  beads  to  monitor  the  condenser 
temperature.  Three  thermoelectric  condensers  were  placed  on 
top  of  the  aluminum  plate.  A  thin  groove  was  cut  around  the 
edges  of  the  inside  face  for  a  0-ring  gasket.  The  plate  was 
screwed  over  the  test  chamber  with  eight  stainless  steel 
screws  of  1/16  inch  diameter  and  compressed  the  0-ring  to  form 
a  watertight  seal.  Figure  (8)  shows  the  cover  from  the  top 
with  the  condensers  on  it. 

B.  INSTRUMENTATION 

Separate  voltage  taps  across  the  heating  element  and 
across  the  precision  resistor  are  used  to  compute  heater 
resistance  and  heating  element  power.  The  heating  element 
temperature  is  computed  from  its  resistance  via  temperature 
calibration  curve  obtained  prior  to  the  experiment. 
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1.  Power  to  the  Heaters 


Each  heater  was  run  in  series  with  a  precision 
resistor  that  was  measured  to  have  a  0.1  ohm  resistance.  All 
four  heaters  were  in  parallel  with  a  Hewlett  Packard  model 
number  6289A,  0-40  volts,  0-1.5  amperes,  DC  power  supply.  Both 
the  voltage  across  the  precision  resistor  and  the  chromel  wire 
were  measured  independently.  The  current  to  each  heater  was 
calculated  by  the  product  of  the  heater  voltage  and  heater 
current.  Both  voltages  were  measured  by  a  Hewlett  Packard 
model  3852S  data  acquisition  system  containing  a  Hewlett 
Packard  model  44701A  digital  voltmeter,  all  controlled  by  a 
Hewlett  Packard  model  300  computer. 

2 .  Temperature  Measurement 

The  thermocouples  for-  the  bulk  fluid,  ambient,  ice 
reference  point  and  coolers  inside  plate  temperatures  were 
referenced  individually  to  an  electronic  Ice  Point  Reference 
as  seen  in  Figure  (9)  .  Each  reference  thermocouple  was 
connected  such  that  its  constantan  lead  was  connected  to  a 
constantan  lead  of  the  measurement  thermocouple.  The  copper 
leads  of  each  measurement  and  reference  thermocouple  were 
connected  directly  to  a  Hewlett  Packard  model  44705A  20 
channel  relay  multiplexer  and  inserted  into  the  data 
acquisition  system.  The  data  acquisition  system  then  measures 
the  ice  referenced  voltage  of  each  thermocouple.  The  voltages 
were  converted  directly  into  temperatures  in  the  controlling 


11 


computer  program  by  using  a  fourth  order  polynomial,  fit  to 
the  thermocouple  manufacturer's  calibration  data  for  copper- 
constantan  thermocouples  (Hazard,  1986) .  The  locations  of  the 
thermocouples  in  the  chamber  are  presented  in  Figure  (10) . 
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Figure  1.  System  Configuration. 


Figure  2 .  Phonograph  of  Tash  J^>paranus . 
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Figure  3.  Isometric  View  of  the  Test  Chamber. 
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Key : 

1 .  Washers 

2 .  Chromel  Wire 

3.  Screw 

4 .  Nut 

5.  Power  Line 


Figure  5.  Sch«matic  Drawing  of  Test  Section  in  Vessel. 
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Figure  7.  Technical  Data  for  Thermoelectric  Coolers. 
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Figure  8.  Cover  Plate  with  Thermoelectric  Coolers. 


Figure  9.  Thermocouple  Connection  Schematic. 


rigur*  10.  Thtt  Locations  of  tha  Tbarmocouplas  in  tha  Chambar . 
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III.  EXPERIMENTAL  PROCEDURE 


A.  NORMAL  OPERATION 

In  carrying  out  the  experiments,  the  procedure  was  as 
follows : 

1.  The  ice  point  reference  device  was  turned  on.  A  mercury 
thermometer  inserted  into  the  ice  base  ensured  that  the  ice 
bath  is  at  0  ±  0.05°C  within  one  hour. 

2.  The  data  acquisition  system  was  then  turned  on.  The 
internal  voltmeter  in  the  data  acquisition  system  requires  one 
hour  warm-up  period. 

3.  The  computer  and  printer  were  turned  on.  The  data 
acquisition  program  contained. in  Appendix  D  was  loaded. 

4.  The  HP  6289A,  0-40  volt,  DC  voltage  supply  was  set  to 
zero  and  energized.  The  current  limiter  on  the  40  volt  power 
supply  was  set  to  the  maximum  current  level. 

A  sampling  of  all  thermocouple  temperatures  was  performed 
and  printed.  The  power  was  increased  incrementally  and  the 
output  from  the  power  calculation  program  was  checked. 

The  bulk  fluid  temperature  as  given  by  the  wire  resistance 
must  match  the  value  given  by  the  bulk  fluid  thermocouples. 
When  the  results  showed  that  the  temperature  stratification  to 
be  less  than  0.1 0°C  the  experiment  was  ready  to  begin.  The 
power  to  the  test  section  is  then  increased  by  controlling  the 
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test  section  heater  voltage.  The  voltage  was  increased  in  0.2 
volt  increments  until  the  initiation  of  boiling  occurred,  at 
which  time  0.1  volt  increments  were  made  until  the  maximum 
attainable  heat  flux  was  reached.  This  occurred  at  5.7  volt 
for  the  power  supply,  making  use  of  the  HP  6289A,  0-40  volt, 
0-1.5  Amperes  DC  voltage  supply  up  to  8.1  volts,  which 
corresponds  to  the  maximum  value.  The  power  supply  has  a  reset 
function  to  protect  the  device  itself.  The  power  is  then 
decreased  in  the  same  fashion.  At  each  power  setting  the 
system  was  permitted  to  stabilize  for  a  few  minutes  and  the 
following  data  were  recorded:  The  heater  voltage,  precision 
resistor  voltage,  two  thermocouple  indications  of 
thermoelectric  cooler  temperatures,  and  three  thermocouple 
readings  of  pool  temperatures.  This  procedure  was  repeated  for 
every  set  up. 


B.  DATA  REDUCTION 

The  test  wire  temperature  was  found  using  the  calibration 
formula  between  thermal  resistance  and  temperature  value  of 
the  chromel  wire. 


fi  =  H„(l+EAr) 

where, 

AT  =  -  Tfc 

Tg  is  the  wall  temperature 

Ti,  is  the  bulk  liquid  temperature 

e  is  the  slope  of  the  calibration  curve. 


24 


Details  of  the  calibration  slope  value  are  presented  in 
Appendix  C.  This  study  utilizes  the  value  of  0.0004  C/ohm 
offered  by  Kaye  (1976) . 

The  starting  value  of  the  wire  resistance  is  accepted  as 
the  Ro  in  the  previous  formula. 

The  heat  flux  is  calculated  from  the  following  formula: 

Flux  =  Power /A 

where,  A  =  n  D  L  is  the  total  wetted  surface  area,  and  D  and 
L  are  diameter  and  length  of  the  wire  respectively. 

The  bulk  temperature  of  the  liquid  was  measured  by  taking 
the  arithmetic  average  of  the  three  thermocouple  readings. 

The  Rayleigh  number  based  on  diameter  is  given  as: 

°  va 

where, 

(p,Cp) 

All  the  properties  in  the  above  equations  are  calculated 
at  the  bulk  temperature. 

Determination  of  Nusselt  number  is: 


where  h 
from 


is  the  heat  transfer  coefficient  and  was  determined 
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^  *  (  T~^ - T 

'  ■‘■s  ■‘■bulk  ' 

where, 

Ts  =  wire  surface  temperature  (°C) 

h  =  heat  transfer  coefficient  (W/m^  •  C) 
g"  =  heat  flux  (W/m^) 

Churchill  and  Chu  (1975)  have  recommended 
correlation  formula  for  the  natural  convection 
horizontal  cylinders. 


Nu^={ 0.60 


0.387  Ray® 
[  1  +  (  0.559/Pr 


It  is  valid  for  a  wide  Rayleigh  number  range 

10“®  <  Ra^  <  10^^ 

where  the  Prandtl  number  (Pr)  is  defined  as 


a  single 
from  long 
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IV.  RESULTS  AMD  DISCUSSION 


Data  were  taken  and  the  related  graphs  can  be  grouped  for 
discussion  purposes  as; 

1.  individually  powered  wires 

2.  pumping  effect  runs 


A.  INDIVIDUALLY  POWERED 

Before  investigating  the  pumping  effect  of  the  bubbles, 
individual  boiling  experiments  on  the  wire  heaters  were 
conducted. 

Figure  (11)  shows  wire  1  (the  lowest  wire)  where  distinct 
single  phase  and  nucleate  boiling  regimes  are  observed.  When 
the  heat  flux  was  increased  in  small  steps,  the  single  phase 
curve  reaches  17°C  superheat.  At  about  20°C  an  abrupt 
temperature  decrease  occurs  as  nucleate  boiling  begins. 

Figures  (11)  and  (12)  are  essentially  the  same.  The 
coordinates  are  heat  flux  versus  surface  temperature  of  the 
chromel  wire.  Figure  (12)  is  on  the  log-log  scale  while  Figure 
(11)  is  semilog  graph. 

Figure  (13)  is  the  single  run  boiling  curve  of  heater 
number  3.  It  was  the  only  one  to  show  overshoot  of  less  that 
7°C. 
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Figure  (14)  for  wire  number  1  shows  the  same  type  of 
behavior  as  Figure  (11)  .  In  this  case  data  were  taken  at 
different  times  in  order  to  demonstrate  the  repeatability  of 
the  experiments.  Data  points  for  two  different  runs  did  not 
coincide,  but  that  is  because  the  bulk  fluid  temperature  is 
not  the  same  for  all  experiments,  but  depends  on  the  ambient 
temperature  of  the  room. 

Figure  (15)  indicates  all  individually  powered  heaters  on 
the  same  graph.  The  difference  in  the  overshoot  may  be  because 
of  the  surface  geometry  of  heater  number  3. 

The  large  aspect  ratio  of  the  heaters  resulted  in  minimal 
conduction  losses  to  the  supports  and  yielded  a  uniform 
temperature  along  the  length. 

The  boiling  curves  for  all  four  wires  are  practically 
superimposable,  even  at  the  departure  from  nucleate  boiling. 
They  differ  only  in  the  heat  flux  at  which  boiling  incipience 
occurred  and  at  which  full  nucleate  boiling  began. 

Figure  (16)  shows  the  best  fit  to  the  natural  convection 
portion  of  the  experimental  data  for  wire  number  1,  and  is 
compared  with  Churchill  (1975)  correlation  values.  There  is  a 
constant  offset  from  the  correlation  values  because  of  the 
wall  efffect  except  for  the  low  Rayleigh  numbers  of  that  run 
which  also  shows  considerable  scatter  in  the  data. 

The  saturation  temperature  of  the  FC-72  is  about  56°C. 
Tsat  “  "^buik  is  the  subcooling.  Figure  (17)  shows  the  effects  of 
subcooling.  At  boiling  incipience,  a  liquid  mass  near  the 


28 


heated  surface  is  slightly  superheated,  even  though  the  bulk 
fluid  is  subcooled.  Each  run  began  with  zero  heating.  The  heat 
flux  value  was  increased  incrementally,  continuing  up  the 
single-phase  natural  convection  line  until  boiling  was 
initiated.  The  figure  indicates  that  as  the  subcooling 
increases  the  curve  shifts  to  the  left,  displaying  significant 
increases  in  the  overshoot  and  heat  fluxes.  With  bulk  liquid 
at  15°C,  it  shows  large  amount  of  superheat,  probably  a  bubble 
begins  to  form  on  the  wire  and  as  soon  as  the  bubble  grows  a 
bit  it  encounters  subcooled  liquid  and  collapses.  The  liquid 
is  so  cold  that  it  doesn't  allow  the  bubble  to  grow.  As  the 
bulk  temperature  of  the  liquid  gets  warmer,  it  allows  the 
bubbles  to  grow.  Subcooling  results  in  smaller  bubbles,  which 
are  less  likely  to  spread  to  neighboring  nucleation  sites  than 
the  larger  bubbles  associated  with  saturated  boiling. 

B.  PUMPING  EFFECT 

For  this  type  of  experiment,  neater  numoer  1  was  used  as 
the  pumper. 

Wire  number  4  was  fixed  at  various  heat  flux  values,  while 
heater  number  1  was  increased  in  a  stepwise  manner. 

Figures  (18)  and  (19)  show  the  pumping  effect  of  the 
bubbles  generated  by  wire  number  1  on  the  temperature  of  wire 
number  4.  Clearly  Figure  (19)  gives  better  understanding  by 
way  of  marked  corresponding  points.  As  the  temperature  comes 
closer  to  the  temperature  at  which  the  boiling  incipience 
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starts,  wire  number  4  is  little  effected.  Just  after  wire 
number  1  is  subjected  to  overshoot  it  also  has  similar  a 
experience,  i.e.,  a  sharp  jump  .  Actually,  it  was  the  general 
pattern  for  most  of  the  other  runs.  This  suggests  that  boiling 
can  be  started  by  the  pumping  effects  of  the  rising  bubbles 
which  trigger  the  heater  located  above  it.  The  characteristics 
of  the  bubbles  leaving  the  lower  heater  are  uncertain.  They 
may  be  the  wake  of  the  buoyancy-induced  motion,  especially 
when  the  upper  wire  was  fixed  at  a  high  heat  value.  A  similar 
trend  can  be  observed  from  Figure  (20)  .  One  thing  to  notice  is 
that  heater  number  4  is  fixed  at  a  heat  flux  that  corresponds 
to  a  starting  value  of  surface  temperature  of  50°C  for  this 
run.  This  is  below  the  temperature  that  meets  the  overshoot. 
Less  jump  in  the  temperature  is  associated  with  the  lower 
fixed  heat  flux  value.  If  we  compare  it  with  Figure  (19), 
which  was  at  higher  heat  value  than  this  one.  Figure  (21) 
provides  a  detailed  look  at  Figure  (20) . 

Figure  (22)  shows  wire  number  3  fixed  at  37,490  W/m^.  It 
has  the  points  marked  in  chronological  order.  The  3rd  wire 
heater  is  located  closer  to  the  first  wire  than  is  the  fourth 
wire,  so  it  allows  most  of  the  bubbles  to  reach  the  upper 
heater  with  collapsing  due  to  surrounding  subcooled  liquid  FC- 
72. 

Figure  (23)  shows  another  combination  of  heater  number  3 
with  heater  number  1.  Heater  number  3  is  kept  at  20,450  W/m^. 
This  was  very  low  heat  flux  value  that  the  triggering  effect 
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did  not  show  itself.  Bubbles  still  caused  some  decrease  in  the 


temperature  of  the  upper  wire. 

Figures  (24a)  and  (24b)  show  very  similar  results.  In  this 
case  wire  number  2  is  )cept  at  20,530  W/m^.  Since  wire  number 
3  is  the  closest  to  wire  number  1.  Being  so  close  helped  it  to 
lose  some  heat  corresponding  to  a  decrease  in  its  temperature. 

Figures  (25a)  and  (25b)  show  no  effect  from  the  bubbles 
from  below.  Here  heater  number  2  is  Icept  at  a  very  high  heat 
flux  of  212,700  W/m^,  which  is  already  boiling.  It  is  not 
sensitive  any  more  to  the  bubbles  produced  by  heater  number  1. 
Consequently  no  jumps  or  apparent  decrease  on  the  surface 
temperature  of  the  number  4  wire  are  seen. 

Figure  (26)  shows  similar  trend  as  the  early  graphs.  In 
this  case  the  bubbles  are  produced  by  wire  number  2  and  the 
effects  on  wire  number  3  are  presented.  Heater  number  3  had 
the  temperature  decrease  as  the  pumper  had  overshoot 
experience . 
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Figure  11.  Wire  Number  1  Individually  Powered,  Semi-log. 
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INDIVIDUALLY  POWERED 


TSURFACE  (C) 

Figure  12.  Niunber  1  Beat  Flux  vs.  Surface  Temperature, 

Log-log. 
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INDIVIDUALLY  POWERED 


■■■■ 

!■■■■ 

i 

■ 

b" 

1 

- B 

B 

"■ 

B 

B 

B 

B 

B 

- B — 

BmHBBmM 

■  ”■  ■ 

■  . . 

LEGEND 
■  =  WIRE  #3 

- 

- 1 

1 - 1 - ! - 1 - 1 - 1 - 1 

23.0  34.0  45.0  56.0  67.0  78.0  89. 

TSURFACE  (C) 


Figure  13.  Number  3  Shows  Less  Orershoot. 
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Figure  14.  Ntunber  1  at  Different  Times. 
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INDIVIDUALLY  POWERED 
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Figure  15.  All  Wiras  Individually  Powarad. 
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Figure  16.  Best  Fit  Xxperisientel  Date  Coiig>»red  with 
Churchill  (1975)  Correlation. 
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Figure  17.  Siibcooling. 
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WIRE  #4  AT  CONS.HEAT  FLUX 
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Figura  18.  Puaqping  Effact  with  Wira  Numbar  4  at  42,800  W/m 


Figure  19.  Pumping  Effect  with  Wire  Number  4 
(Detailed  Figure) . 
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WIRE  #4  AT  CONST.  HEAT  FLUX 
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Figura  20.  Nira  Numbar  4  at.  Constant  Baat  Flux  (38,600  W/m^)  . 


41 


WIRE  1  &  4(FIXED) 


rigur«  21.  Detail  of  Flgura  20. 
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HEAT  FLUX  (W/M  2) 


WIRE  1  &  3  (FIXED) 


rigurs  22.  Wlra  Nuiiib«r  3  at  Constant  Boat  Flux  Valua  of 

37,490  1f/m^ 


43 


HEAT  FLUX  (W/M  2) 
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WIRE  #3  AT  CONS.  HEAT  FLUX 
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WIRE  #2  AT  CONS.  HEAT  FLUX 
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Figure  24a.  Wire  Mumber  2  at  20,530  W/m^ 


45 


WIRE  #2  AT  CONS.  HEAT  FLUX 
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Figure  24b.  Detail  ot  Figure  24a. 


46 


WIRE  #2  AT  CONS.HEAT  FLUX 
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Flguro  25a.  Wire  Number  2  Kept  at  212,700  W/m^. 
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WIRE  #2  AT  CONS.HEAr  FLUX 


rigur#  25b.  Detail  of  Figure  25a. 
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HEAT  FLUX  (W/M  2) 

4 
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Pigur*  26-  Wire  Number  2  with  Wir«  Number  3  »t 
Constant  Boat  Plux. 
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V.  CONCLUSIOHS 


The  figures  give  the  idea  of  the  pumping  effects.  When  the 
pumping  heater  is  powered  progressively,  the  upper  heater  will 
react  as  indicated  in  the  following. 

•  If  the  upper  wire  is  located  close  to  the  lower 
numper,  it  will  meet  most  of  the  bubbles.  It  will  lose  some 
heat  depending  on  how  close  it  is. 

•  The  initial  heat  flux  value  of  the  upper  heater  is 
also  important.  The  amount  of  temperature  jump  which 
simultaneously  occurs  when  the  pumper  has  overshoot  depends  on 
how  warm  the  upper  wire  is. 

•  If  the  upper  heater  already  has  bubbles,  this  means  it 
had  incipience  before,  and  the  bubbles  will  not  help  it  lose 
more  energy. 
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VI.  RECOMMEMDATIONS 


1 .  Use  platinum  wires  rather  than  chromel  wires  to  get 
better  temperature  measurement. 

2.  Before  the  start  of  each  run,  boil  the  pool  for  a 
sufficient  time  to  degas  the  liquid. 

3.  Modify  the  software  for  taking  data  by  adding  a  WAIT 
command  to  obtain  the  time  function. 

4.  An  addition  of  a  foil  type  heater  placed  horizontally 
at  the  bottom  of  the  chamber  may  help  create  a  uniform 
temperature  throughout  the  bulk  fluid. 

5.  An  extension  of  tygon  tubing  of  about  one  foot  length 
would  take  care  of  overflow  and  hydrostatic  head  determined  by 
the  liquid  level  inside  the  filling  tank. 

6.  The  present  study  deals  with  only  subcooling  effects. 
Saturated  pool  boiling  should  also  be  studied.  In  order  to 
have  constant  temperature  in  the  bulk  fluid,  the  thermo¬ 
electric  condensers  and  the  bulk  heaters  can  be  run 
simultaneously . 
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APPENDIX  A.  SAMPLE  CALCOLATIONS 


A.  DETERMINATION  OF  INPUT  POWER 

For  the  sample  calculations,  the  output  values  of  wire 
number  1  were  used: 

Power  ( 1 )  =  ( 11 ) 

where  volt  (12)  and  volt  (11)  are  the  voltages  across  the 
precision  resistor  and  wire  respectively. 

=  I  j(1.988)  =  0.96  watts 

B.  CONVERSION  FROM  VOLTAGES  TO  TEMPERATURES  ("C)  FOR 
THERMOCOUPLES 

The  HP  3852  A  Data  Acquisition  System  channels  1-9  record 
the  voltage  readings  in  millivolts.  The  following  polynomial 
formula  is  for  copper  constantan  thermocouples. 

T{°C)  =  0.10086091  +  (  25727.9  ♦  emf)  -  (  767343.5  •  emf^ ) 

+  (  7802.5596  ♦  emf^  )  -  (  9247486589.6  ♦  emf‘^)  + 

(  6.98F1I  ♦  emf^  )  -  (  2.66£13  *  emf^ )  +  (  3.94FI4  *  emf^  ) 

Using  thermocouple  1 

emf  =  0.000779  volt  r=19.72°C 

C.  CALCULATION  OF  WIRE  TEMPERATURE 

The  calibration  curve  for  the  chromel  wire  (Kaye,  1976) 
yields  the  slope  e  of  4.10"^  ohm/°C. 
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=  [  1  +  e  r- rj 

Here  and  R  are  the  resistances  of  the  wire  at  the  reference 
point  and  at  any  point  to  be  measured  Rg  and  ohm(l)  are  the 
corresponding  variables  used  in  the  software. 

Ts(l)  =  (((ohm(l)/RJ  -  1)  /  Kaye  )  + 

where, 

Ts(l)  =  wire  temperature 

Kaye  =  0.0004  ohm/°C  constant  slope 

To  =  wire  reference  temperature  (°C)  . 

TJl)  =  (( (4.1009/4.075)  -  1)  /  4E-4)  +  19.64 

=  35.5°C 

D.  CAIiCniATION  OF  THE  AVERAGE  BULK  TEMPERATURE 

The  average  temperature  of  the  liquid  was  simply  the 
arithmetic  average  of  the  three  thermocouple  values  immersed 
in  the  bulk  fluid. 

Ttuxk  =  (T(l)  +  T(2)  +  T(3))  /  3 

=  (19.9  +  19.54  +  19.6)  /  3  =  19.68°C 


DETERMIMATIOM  OF  HEAT  FLUX 


where, 


HFLUX(l) 


AREA 
D,  L 
AREA 


TC  *  D  *  L  =  (3.1415926)  D  L 
Diameter  and  length  of  the  wire. 
(3.1415926)  (0.005)  (0.0254^)  (3) 


=  3.04E-5 


HFLUX-  =77460 

3.04^-^ 
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r.  FLUID  PROPERTIES  DETEBMIHXTIOM 

1.  Thermal  Conductivity  k(1f/iB*‘C) 

From  Figure  5  of  the  3M  Corporation  Flourinert  Product 
Manual,  the  thermal  conductivity  coefficient  curves  have  been 
determined  to  be 

JC=  (  0.6033  -  0.00115  *  r^)  /  10  .^^  =  0.lJ^ 

at  r*  =  50°C  it  yields  0.0545  W/m. 

2.  Liquid  Density  p  (kg/m^) 

Using  the  expression  in  Table  4B  and  constants 
presented  in  Table  4C  of  the  Product  Manual  yields 
()=  (1.740  -  0.00261  TJ  *  1000 
Tb  film  temperature  must  be  in  units  of  °C.  At  T^  =  50°C  this 
yields  p=  1609.5  kg/m^. 

3.  Kinematic  Viscosity  v  (mVs) 

From  Figure  (3)  and  determining  a  4th  order  curve  fit 

yields 

V  =  (0.54  -  0.0116  ♦  Tfc  +  0.0002085  *  T/;  •  lE-6 
at  Tt,  =  50°C  this  yields  v  =  4.8E-7  mVs. 

4.  Density  of  Vapor  p^  (kg/m^) 

From  Table  13A  of  the  Product  Manual,  values  for  vapor 
density  can  be  taken. 

For  the  range  of  20-80°C  the  data  can  be  approximated 
as  a  fourth  order  polynomial. 

Vv  =  (246  -  10.93  *  Tt,  -h  0.213  *  tj  - 

0.00193  *  Tt,^  +  6.55E-6  *  T^,*)  /  338 
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Tt  =  50°C  it  yields 

Py  =  10.67  (kg/m^ 


where 


5.  Specific  Heat  Co  (J/kg"C) 

For  all  Flourinert  Electrochemicals  (Figure  4) 

Cp  =  (0.241111  +  3.70374E-4  *  Ty,)  *  4186 
at  T  =  50°C.  Calculation  yields  1086.7  J/kg°C. 

6.  Thermal  Diffusivity  a  (m‘/s) 

a  =  k/pCp 

at  T  =  50°C/  calculation  yields  3.115E-8  m^/s. 

7.  Thermal  Expansion  Coefficient  P  (1/‘’C) 

Table  4B  in  the  Product  Manual  shows 

p  =  0.00261  /  (1.740  -  0.00261)  *  Ty,) 
at  r  =  50°C.  Calculation  results  in  0.00162  (1/°C. 

8.  Prandtl  Kumber  (Pr) 

Pr  =  v/a 

at  T=50°C.This  yields  Pr  =  15.4 

9.  Calculation  of  Latent  Heat  hfg  (J/kg) 

Table  13A  of  the  Product  Manual  shows  the  values  can 
be  approximated  by  a  linear  equation. 

hfg  =  (8057  -  22.73  *  Ty,)  *  10.79 
where  10.79  is  the  conversion  factor  from  Cal/Mole  to  (J/kg) . 
At  T  =  50°C  calculation  yields  74672  (J/kg) . 
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APPENDIX  B.  ONCBRTAINTY  ANALYSIS 


A.  UNCERTAINTY  IN  AREA 

(Oc  =  .0001  inches 

cOi  =  ±0.010  inches 

A  =  n  D  L 


where : 
D 
L 
A 


Wire  diameter 
Length  of  heater  wire 
(3.141593  •  0.005  *  3) 
0.0471  in^ 

3.04  *  10'^ 


Ba 

-373 


(n  L) 


Ba 


(n  D) 


Therefore, 

(0^  =  [  <  ^  ^ 
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1  0.0001  .  / 

0.020  \0 

H  0.005  /  \ 

“3-/  J 

=  0.021  or  2.1% 

A 

0)  =  6.38  X  10-*  in^  =  9.89  x  lO”" 

B.  UNCERTAINTY  IN  POWER 

=  ±0.05  Ampere  Taken  from  (Hazard,  1986) 

calculation) 

=  ±0.5  volt 

=  VI 

=  I 

=  V 


<Oi 

OK 

Q 

dQ 

do 


cob  =  [  ( 


=  0.059  or  5.9% 
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C.  OHCKRTAIIlTy  IN  SEAT  FLUX 


<? 


Q 

A 


dq  _  1 

"TS  A 


dq 

■93 


(Oq 


cog 

<7 


/■  \2 

/  \2  1 

“a 

+ 

\  ^  J 

=  [  0.021^  +  0.059^ 


0.062  or  6.2% 


D.  UNCSRTAINTY  IN  TEMPERATURE 


Using  thermocouples: 
Variable 

Voltmeter  resolution 

Ice  reference  bath 
Temperature 

Polynomial  Temperature 

Conversion 

R 


Uncertainty 

0.026°C 

1.0  mv  resolution 
0.05°C 

0.00663°C  RMS 

1.0% 


=  [  (0.0025)^  +  (0.05)^  +  0.00663 

(OTt,  =  0.063  or  6.3% 


E.  UNCERTAINTY  OF  WIRE  SURFACE  TEMPERATURE 

From  the  calibration  formula 


where, 

AT 

e 

Tt 

From  (1) 


R  =  R,,  [  1  +  e  AT  ] 


Ts  - 

Slope  of  the  calibration  curve 
Bulk  fluid  temperature. 


(1) 


T 


S 


/  e 


For:  Ts  =  44.14°C 
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60 


-1 


^  "• 


9.43 
44.14  ) 


.063  \- 

Tr^l 


0.21  or  21% 

0.21  (44.14  -  23.99) 


=  4 . 3°C 


6.  UNCERTAINTY  IN  HEAT  TRANSFER  COEFFICIENT 


q  is  the  heat  flux 


q  =  20190  W/m^ 
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dh 


Af’ 


AT  =  44.14  -  23.99 


<A, 


=  20.15°C 


©h 

“?r 


CO. 


V 


_q 


/'  V 

0)at 

~Sf 


1/2 


(  0.062  ) 


=  0.222  or  22% 
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APPENDIX  C.  CALIBRATION  OF  CHROMEL  WIRES 


The  chromel  wires  of  about  3  inch  length  were  immersed  in 
a  Rosemount  constant  temperature  distilled  water  bath,  and 
were  compared  with  a  platinum  resistance  thermometer  as  a 
standard.  The  temperature  was  raised  from  30°C  to  90°C  in 
increments  of  about  five  degrees-  At  each  increment,  the 
temperature  was  allowed  to  stabilize  for  five  minutes,  then 
measured  using  a  Newport  digital  thermistor  devise  for  the 
0.005  inch  diameter  chromel  wire.  All  recorded  data  are 
presented  in  the  following  table.  The  corresponding  best  fit 
line  is  also  presented  on  the  graph.  Two  sets  of  data  were 
taken  at  different  times. 


Thermistor 

(C) 

T/C 

Resistance 

(Ohm) 

29 

30.2 

4.74 

38.3 

39.7 

4.76 

49.5 

51 

4.78 

61.9 

63.4 

4.8 

76.2 

77.2 

4.83 

87.7 

88.4 

4.85 

52.3 

53.8 

4.79 

57.5 

59.1 

4.8 

60.5 
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4.8 

68.3 

69.6 

4.81 

70.7 

72 

4.82 

71.3 

72.4 

4.82 

86.2 

86.9 

4.84 
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JWIPL 
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4.70 


4.72  4.74  4.76  4.78  4.80  4.82  4.84  4.86 

Resistance  (Ohm) 


Figur*  Cl.  ChrooMl  Wir*  Calibration. 


APPENDIX  D.  SOFTNARE 


Software  written  in  the  BASIC  language  converts  voltage 
measurements  into  temperature  values  for  the  bulk  fluid,  ice 
reference,  and  coolers  inside  the  thermocouples.  It  uses  the 
a  value  of  0.0004  Ohm/°C  to  calculate  the  surface  temperature 
offered  by  Kaye  (1976) .  The  printed  output  covers  bulk  fluid, 
wire  surface  temperatures  and  heat  flux. 

A  FORTRAN  program  reads  the  output  values  created  by  the 
previous  program  and  calculates  liquid  properties  R^,  Nu,  and 

The  GRAPHER  program  graphs  the  data  and  prints  it  using  a 
HP-7440A  Plotter. 

Initial  experiment  data  was  fed  into  EASYPLOT  software  on 
the  Naval  Postgraduate  School's  main  frame  computer  via  floppy 
diskettes . 
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n  OPTION  BASE  1 
\2  BEEP 

t3  PRINT  ‘INPUT  DATE* 

U  INPUT  TarihS 

15  PRINT  •  INPUT  TINE  ‘ 

16  INPUT  Zan 

IB  IPRINT  ‘SCANNIN6  NEA5* 

19  REAL  Enf(7),T(7>.0<8).Uolt( 11 tl9).Power<4).Hflux(4)jB(4>,0hn(i) 
2B  REAL  Suno 

25  OUTPUT  709i'CONFHEAS  OQU ,101 -107 .USE  0* 

26  FOR  I-l  TO  7 

27  ENTER  709iEmf(n 

29  IPRINT  ItEinf<l) 

30  NEXT  I 

31  OUTPUT  70S«’CONFHEAS  DCV . 1 11 -1 1 9 .USE  0* 

32  FOR  I-Il  TO  19 

33  ENTER  709tUolt<I) 

35  IPRINT  IiUoltd) 

3B  NEXT  I 

37  0< I  )-. 10096091 

36  0(2 >-25727.9 
3S  0<  3  >—767345.9 

40  0(4 >-76025595 

41  0(5  >—9247466589 

42  0(5  >-6.95E-*'n 

43  0(7)— 2.665+13 

44  0(6>-3.94E+14 

45  FOR  I-l  TO  7 

46  Suno-0 . 

47  FOR  J-I  TO  8 

48  Suno-3uno+0( J  )»Enf( I >"( J-1  > 

49  NEXT  J 

50  T(I>-Suno 

51  Tb-(T(2>+T(3>+T(l  >>/3 

52  NEXT  I 

53  PRINT  'DATE  j'jTarih* 

54  PRINT  'TIME  »  ";Zan 

55  PRINT  'T'd.Td);'  C  UPPER  PORTION  OF  THE  FLUID’ 

56  IPRINT  USING  'T' t I .DD.D(T( I >i *  C  UPPER  PORTION  OF  THE  FLUID" 

57  PRINT  ■T’«2,T(2)«'  C  LOWER  PORTION  OF  THE  FLUID’ 

59  PRINT  ‘T'jS.TISli’  C  LOWER  PORTION  OF  THE  FLUID' 

59  PRINT  •T"j4.T(4)i‘  C  TE  COOLER  INSIDE  AL.  PLATE’ 

60  PRINT  *T";5,T(5)r  C  TE  COOLER  INSIDE  AL.  PLATE" 

61  PRINT  ‘T'jS.TISIj*  C  AMBIENT  (ROOM  TEMPERATURE)' 

TEMPERATURE  )• 


52  PRINT  *T'}7,T(7>;“C  ICE  (REFERENCE  TEMPERATURE >■ 

53  (NEXT  I 

54  PRINT  •TBULK-";Tb 

55  Ipouer  coiculatlons»**»»***»»»» 

70  PoMerd  >-<Uolt(12)/.l)»Uolt(Tri 

50  PoMer(3)-<Volt<I7)/.l 15  > 

51  theater  surface  area  calculations 

52  !  ••«•*•••••••••••••••••••••••••••• 

53  Dla> .005* .8254 

94  Al-3>.e254 

35  Alan»<3. 1415325 >*Oia*Al 
37  Iheat  flux  calc. 

93  . . 

99  Hflux< I )“Pouer< 1 >/Alan 

100  Kf lux< 2 )*Pawer< 2 )/Alan 

101  Hflux(3 )"Pouer( 3  i/Alan 

102  PRINT  "power  1 iPouer< 1  i 

103  IFRINT  ’pouer2“’ ;Pouer(2 1 

104  'PRINT  *power3“" ;Pouer(3 ) 

105  I . * . . 

107  Ohfi(  1  j-Uoltdl  >/(Uoit(i2>/.i  > 

103  !0hnv2)-Uoltd4)/(Uoltd6)/.l  ) 

109  0hn(3)-UoltdB)/(Uoit(  17)/.1  ) 

111  tPRINT  •0HH2-*i0hn<2) 

112  tPRINT  "0HMI-"  iOhi<i(  1  ) 

114  IPRINT  "OHMi-’jOhnlj) 

1 1 3  0t*Tb-7a 

115  !  Ta  IS  THE  TEMP  OF  BULK  FLUID  Tb  WHEN  POWER  OFF 
117  !  Ro  IS  THE  RESIS.  Or  CHROMEL  WIRE  ACTUALLY  IT  IS 

113  i  THE  Ohnd  )  VALUE  WHEN  THE  BULK  FLUID  TEMP  (Tb-Tamb)  AND  POWER  OFF 

119  1  <  POWER  OFF  MEANS  Pouer( 1 )"0.0084  and  it  is  obtained  by  setting  ) 

120  !  POWER V «)* 0.0038 

121  !  R0I-4.0S3  R04-3.7981 

123  !  (VIZ  40  V  Supplyist  to  00.1  0.10  position  WHEN  1&4  INSERTED 

124  !  Kaye  is  the  value  of  the  calibration  slope  ReflS.M  YOU  et  all 

125  I  . . . . 

125  Ta*Tb 

127  Ro1>4.05S 

128  Ro3"3.7381 

129  Kayeo.0004 

130  R"4.053»( 1  +  .0004*Ot  1 

131  Ts(  1  )"( ( <  Ohnv  1  )/Ro  1  >-l  )/Kaye  )+7e 

132  Ts(  3  )"( ( (  Ohn(  3  )/Ro3  )-1  i/Kaye  )+Ta 

133  Tf iln«( 7b+7s( 4  )  )/2 

134  PRINT  "TSl-’iTsI  I  ).’  HFLUX  j'jHfluxd) 

135  PRINT  "TSS-" ;Ts(3)  .*  HFLUX3  »'iHflux(3> 

135  IflRINT  ’T  flln  ;*!7filn 

340  END 
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PARAMETER  (  N>E8  ) 

INTEGER  I 

REAL  TS  (N>  ,  TB  (N)  ,  HF  (N)  ,  NU  (N)  ,  RA  (N>  ,  NUC  (N) ,  TF  (N)  ,  TFI  (N)  ,  BETA  (N) 
,  DENF  (N)  ,  DENV  (N>  ,  KINW  (N)  ,  CP  (N> ,  K  (N)  ,  HFQ  <N>  ,  ALF  (N)  ,  PR<N) 

,  SURT  <N) ,  RAYD  <Nl ,  RAYFB  <NI ,  NUSEX  (N)  ,  NUSC  <N)  ,  OIA,  S,  V  (M) 

,  MU  (N>  ,  Y  (N)  ,  X  (N)  ,  US  (N)  ,  TSAT,  CSF  (N> 


*  READ  T*urfae«,H««t  Flux,Tbulk  FROM  DATA  FILE  • 

*  WRITE  Tsurfaca, Tbulk, H«at  Flux, Nussalt, Raylaigh, Nussalt  (Chur)* 

*  Tfila  ,  bata, density, viscosity, Cp, conductivity, latant  haat, 

*  Prantl, surfaca  tansion 


OPEN  (UNIT»1 ,  FILE*' INPUT. DAT* , STATUS-'  OLD* ) 

OPEN  <UNIT-2, FILE*' OUTPUT. DAT'  ,  STATUS*'  NEW* ) 

OPEN  <UNIT»3, FILE**  PROPEl. DAT* , STATUS**  NEW* ) 

OPEN  (UNIT-4,  FILE**  PROPEE.  OAT*  ,  STATUS*' NEW* ) 

OPEN  (UNIT-5,  FILE**  ROHSE. DAT* , STATUS-'  NEW* ) 

DO  IS  I  -  1,N 

READ  (1,*)  TS(I),HF(I),TB(I) 

TF(I)-(TS(I)*TB(1) ) /a. 

TFl (I)-TF<I)+a73. 15 

BETA  ( I )  *i3.  0Sa6  i  /  ( 1 .  740-8.  00a6 1*78  ( I )  > 

PRINT*, 'BETA  IS  *',BETA(I) 

DlA-0. 005*0. 0a54 

G«^.ai 

TSAT-56. 

US ( I ) -TS ( I ) -TSAT 

DENF  ( I ) « ( 1 .  740-0. 00a& 1 *T8 ( I ) >  *1000. 

y  ( I )  * (346- 10.  93*TB< I ) *0. ai3*TB ( I ) **a-l9. 3E-4*TS ( 1 ) **3+ 
*  6.55E-6*T8(I)**4)/33a. 

DENV(I)-1. /V(I) 

KINV ( 1 ) - (0. 54-0. 01 16*T8 (I)-Md. 0003005*78 ( I ) **S) *lE-6 
CP  ( I )  -  <0.  341 1  ll-fS.  7037E-4*TB(  I)  >  *4106. 

K < I ) - (-0. 001 15*TF ( I ) +0. 6033) / 10. 

HFG ( I ) » <  a057-33. 73*T8 ( I ) ) *10. 70 

ALF  < I ) *K ( I ) / (DENF ( I ) *CP ( I ) ) 

PR ( I ) -K INV  < I ) /ALF  < I ) 

SURT(I)*(-0, 13*TF<I)+15)/1000. 

MU ( I ) -KINW ( I ) *DENW ( I ) 


C 


c 

c 


30 

30 

40 

50 

15 


Y  ( I )  *  ( HF  <  I )  /  ( MU  ( I )  *HFG  ( I )  ) )  *SaRT  ( SURT  ( I  >  /  (G*  (DENF  ( I )  -DENV  ( I )  )  )  ) 
X  ( I )  *CP  ( I )  *  (TS  ( I )  -TB  ( I ) )  /  (HFG  ( I )  *PR  ( I )  **1.  7) 

CSF (I)-X(I)/(Y(I)  **0. 333) 


CALCULATE  RAYLEIGH  and  NUSSELT  NUMBERS 

RAYD(I)*Q*BETA(I)*(TS(I)-TB(I) )*DIA**3/(KINW(I)*ALF(I) ) 
RAYFB  ( I ) -Q*BETA ( I ) *HF ( I ) *DIA**4/ (K ( I ) *K INV ( I ) *ALF ( I ) ) 
NUSEX(1)-HF(I)*DIA/(K(I)*(TS(I)-TB(I)  )  ) 

NUSC  ( I ) «  (O.  60*  (0.  3a7*RAY0  ( I)  ** ( 1.  /6.  )  /  (  ( 1+  (0.  559/PR(  I ) )  **  ( 
9. /16. ) )**(a. /a7. ) ) ) )**3 


WRITE  (3,30)  TS(I),TB(I),HF(I),NUSEX(I),RAYD(I),NUSC(I) 
WRITE  (3,30)  TF(I),BETA(I),DENF(I),DENW(I),KINW(I),CP(I) 
WRITE  (4,40)  R(I),HFG(I),ALF(I),PR(I),SURT(I) 

WRITE  (5,50)  Y(I),X(I),CSF(I) 

FORMAT  (lX,Fa.  3,3X,F9.1,3X,F13.3,ax,Fa.  3,aX,F13.7,  3X,F15.3  ) 
FORMAT  (lX,Fa.  4,  3X,F9.4,3X,F9.  l,3X,Fa.  4,3X,E10.3,3X,F15.3  ) 
FORMAT  (lX,Fa.  4,  3X,Fl4.4,3X,E10.3,aX,F14.  4,aX,F9.5  ) 

FORMAT  (lX,Fa.  4,  3X,F14.  4,3X,E10.3  ) 

CONTINUE 


68 


LIST  OF  REFERENCES 


3M  Corp.,  Flourinert  Product  Manual,  Conunercial  Chemicals 
Division:  St.  Paul,  Minnesota,  1987. 

Bar-Cohen,  A.,  and  T.  Simon,  "Wall  Superheat  Excursions  in  the 
Boiling  Incipience  of  Dielectric  Fluids, "  Heat  Transfer 
Engineering,  Vol.  9,  no.  3.,  pp.  19-31. 

Bar-Cohen,  A.,  and  T.  Simon,  Experiments  on  Boiling  Incipience 
with  a  Highly-Wetting  Dielectric  Fluid:  Effects  of  Pressure, 
Subcooling,  Dissolved  Gas  Content. 

Benedict,  T.,  An  Advanced  Study  of  Natural  Convection 
Immersion  Cooling  of  a  3x3  Array  of  Simulated  Components  in  an 
Enclosure  Filled  with  Dielectric  Liquid,  Master's  Thesis, 
Naval  Postgraduate  School:  Monterey,  California,  June  1988. 

Berenson,  P.J.,  "Experiments  on  Pool-Boiling  Heat  Transfer," 
International  Journal  of  Heat  Mass  Transfer,  Vol.  5,  pp.  985- 
999,  1962. 

Bergles,  A.,  N.  Bakhru,  and  J.  Shires,  "Cooling  of  High  Power 
Density  Computer  Components,"  Department  of  Mechanical 
Engineering  Projects  Laboratory  Report,  Massachusetts 
Institute  of  Technology,  1968. 

Bergles,  A.  and  Chyu,  "Nucleate  Boiling  from  Porous  Metal 
Coatings,"  Energy  Conservation  via  Heat  Transfer  Enhancement, 
Department  of  Energy  Publication  COO-4649-10,  pp.  5-13,  1979. 

Chowdhury,  S.K.R.,  and  R.H.S.  Winterton,  "Surface  Effects  in 
Pool  Boiling, "  International  Journal  of  Heat  Mass  Transfer, 
Vol.  28,  pp.  1881-1889,  1985. 

Churchill,  S.W.,  and  H.H.S.  Chu,  "Correlating  Equations  for 
Laminar  and  Turbulent  Free  Convection  from  a  Horizontal 
Cylinder,"  International  Journal  of  Heat  Mass  Transfer,,  Vol. 
18,  p.  1049,  1975. 

Frost,  W.,  and  G.S.  Dzakowic,  "An  Extension  of  the  Method  of 
Predicting  Incipient  Boiling  on  Commercially  Finished 
Surfaces,"  ASME  paper  67-HT-61,  1967. 

Hazard,  J.,  Single  Phase  Liquid  Immersion  Cooling  of  Discrete 
Heat  Sources  on  a  Vertical  Channel,  Master's  Thesis,  Naval 
Postgraduate  School:  Monterey,  California,  December  1986. 


69 


Incropera,  F.P.,  and  D.P.  DeWitt,  Introduction  to  Heat 
Transfer^  (John  Wiley  &  Sons,  1976) . 

Kaye,  G.W.C.  and  T.H.  Laby,  Tables  of  Physical  and  Chemical 
Constants  and  Some  Mathematical  FunctionSr  14th  ed.,  (London: 
Longman),  p.  104. 

Lee,  T.Y.,  and  T.W.  Simon,  "High-Heat-Flux  Forced  Convection 
Boiling  from  Small  Region, "  Heat  Transfer  in  Electronics,  ASME 
HTD-Vol.  Ill,  pp.  7-16,  1989. 

Lepere,  V.,  Nucleate  Pool  Boiling  of  High  Dielectric  Fluids 
from  Enhanced  Surfaces,  Master's  Thesis,  Naval  Postgraduate 
School:  Monterey,  California,  December  1980. 

McAdams,  W.H..,  W.E.  Kennel,  C.S.  Minden,  Rudolf  Carl,  P.M. 
Picornell,  and  J.E.  Dew,  "Heat  Transfer  at  High  Rates  to  Water 
with  Surface  Boiling, "  Industrial  and  Engineering  Chemistry, 
vol.  41,  No.  9,  1949. 

Mudawar,  I.,  and  D.E.  Maddox,  Critical  Heat  Flux  in  Subcooled 
Flow  Boiling  of  Fluorocarbon  Liquid  on  a  Simulated  Electronic 
Chip  in  a  Vertical  Rectangular  Channel. 

Nishikawa,  K.,  and  T.  Ito,  "Augmentation  of  Nucleate  Boiling 
Heat  Transfer  by  Prepared  Surfaces,"  Japan-United  States  Heat 
Transfer  Joint  Seminar,  Tokyo,  Japan,  1980. 

Okuyama,  Kunito,  and  Yoshihiro  Ida,  Transient  Boiling  Heat 
Transfer  Characteristics  of  Nitrogen. 

Rohsenow,  W.M.,  "Boiling,"  Handbook  of  Heat  Transfers,  W.M. 
Rohsenow  and  J.P.  Hartnett,  eds.,  (New  York:  McGraw  Hill, 
1973),  chapter  13. 

Tang,  W.,  A.  Bar-Cohen  and  T.  Simon,  Thermal  Transport 
Mechanisms  in  Nucleate  Pool  Boiling  of  Highly-Wetting  Liquids. 

Torikai,  K.,  H.  Shimamune,  and  T.  Fujishiro,  "The  Effects  of 
Dissolved  Gas  Content  upon  Incipient  Boiling  Superheats," 
Proceedings  of  the  4th  International  Heat  Transfer  Conference, 
Vol.  V.,  B2.11,  1970. 

Yilmaz,  S.,  J.  Hwalek,  and  J.  Westwater,  "Pool  Boiling  Heat 
Transfer  Performance  for  Commercial  Enhanced  Tube  Surfaces, " 
ASME  Paper  No.  80-HT-41,  National  Heat  Transfer  Conference, 
Orlando,  Florida,  1980. 

You,  S.M.,  "Boiling  Heat  Transfer  with  Highly-Wetting 
Dielectric  Fluids,"  Ph.D.  Thesis,  University  of  Minnesota, 
1990. 


70 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1 .  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Mechanical  Engineering  Curricular  Office,  Code  ME  1 

Naval  Postgraduate  School 

Monterey,  California  30340 

4.  Professor  M.D.  Kelleher,  Code  ME/Kk  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5000 

5.  Professor  Y.  Joshi,  Code  ME/Ji  1 

Naval  Postgraduate  School 

Monterey,  California  93943-5000 

6.  Naval  Weapons  Suppport  Center  1 

Code  6042 

Crane,  IN  47522 

7.  Dz  K.  Komutanligi  1 

Okullar  ve  Kurslar  Dairesi 

Bakanliklar,  Ankara,  Turkey 

8.  Deniz  Harp  Okulu  1 

Kutuphanesi 

Tuzla,  Istanbul,  Turkey 

9.  Istanbul  Teknik  Oniversitesi  1 

Makine  Fakiiltesi  Kytyphanesi 

Istanbul,  Turkey 

10.  LTJG  Ali  Sukrii  Eren  2 

Harzem  Sok  5/3,  80650 

Celiktepe,  Istanbul,  Turkey 

11.  LT  R.A.  Egger  1 

SMC  #1039 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 


72 


